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The paper extends the result of Harman and Pronzato [Stat. & Prob. Lett., 77:90-94, 
2007], which corresponds to p = 0, to ah strictly concave criteria in Kiefer's (pp-class. Let ^ 
be any design on a compact set ^ C M™ with a nonsingular information matrix M(^), and 
let 6 be the maximum of the directional derivative Fci,^{£,,x.) over all x e We show that 
any support point x, of a (/)p-optimal design satisfies the inequality i^0p(^,x,) > hp[M.{^),d], 
where the bound /ip[M(^), S] is easily computed: it requires the determination of the unique 
root of a simple univariate equation (polynomial when p is integer) in a given interval. The 
construction can be used to accelerate algorithms for (/)p-optimal design and is illustrated on 
an example with A-optimal design. 
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1 Introduction and motivation 

For ^ a compact subset of M™, denote by H the set of design measures (i.e., probability 
measures) on ^ and by M(^) the information matrix 



We suppose that there exists a nonsingular design on ^ (i.e., there exists a ^ G H such that 
M(^) is nonsingular) and we denote by the set of such designs. We consider an optimal 
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design problem on J?r defined by t he maximiza t ion of a design criterion 0(^) = <I>[M(^)] with 
respect to ^ G H. One may refer to iPukelsheimI () 19931 . Chap. 5) for a presentation of desirable 
properties that make a criterion $(•) appropriate to measure the information provided by ^. 
He re we shall fo cuss our attention on design criteria that correspond to the (/>p-class considered 
by iKieferl (jl974l ). More precisely, we consider the positively homogeneous form of such criteria 
and, for any M G M, the set of symmetric non-negative definite m x m matrices, we denote 

-i/p 



$;(M) 



1 



m 



tr(M-P) 



il) 



with the continuous extension $p (M) = when M is singular and p > 0. A design measure 
that maximizes (l)p{(,) = $p [M(,^)] will be said 0p-optimal. Note that when p ^ the 
maximization of $jJ'(M) is equivalent to the minimization of [tr(M~^)]^^''', and thus to the 
minimization of tr(M~*') when p is positive. A classical example is A-optimal design, which 
corresponds to p = 1. Taking the limit of $p (•) when p tends to zero, we obtain <I>^(M) = 
[det(M)]"'^/™' , which corresponds to D-optimal design. The limit when p tends to infinity gives 
<I>oo(M) = Amin(M), the minimum eigenvalue of M, and corresponds to E'-optimal design. Some 
basic properties of i?i>p-optimal designs are briefly recalled in Sect. [2l 

Classical algorithms for optimal design usually apply to situations where ^ is a. finite set. 
The performance of the algorithm (in particular, its execution time for a given required precision 
on (/>(•)) then heavily depends on the number k o f elements in ^ . The case of -D-optimal design 
has re t ained muc h at tent i on, se e, for instance. lAhipasaoglu et al. ( 2008 ) 



Harman and Pronzato ( 



Todd and Yildirim 



20071 ) show how any nonsingular de- 



(120071), |Yu| (j2O10|) andlYul ((20111). 

sign on ^ yields a simple inequality that must be satisfied by the support points of a L)-optimal 
design ^q. Whatever the iterative method used for the construction of ^g, this delimitation of 
the support of permits to reduce the cardinality of ^ along the iterations, with the inequality 
becoming more stringent when approaching the optimum, hence producing a significant accel- 
eration of the algorithm. Put in other words, the delimitation of the support of an optimal 
design facilitates the optimization by focussing the search on the useful part of the desig n space 
The objective of the paper is to extend the results in iHarman and Pronzato! (|2007l ) to the 
0p-class ([1]) of design criteria. The condition obtained does not tell what the optimum support 
is, but indicates where it cannot be. 

The paper is organized as follows. Section [2] recalls the main properties of (pp-optimal design 
that are useful for the rest of the paper. The main result is derived in Sect. [3] and illustrative 
examples are given in Sect. HI Finally, Sect. [5] concludes and indicates some possible extensions. 
The technical parts of the proofs are given in appendix. 
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2 Some basic properties of 0p-optimal designs 



The criteria (•) defined by ([T]) satisfy ^p{Im) = 1 for Im the m-dimensional identity matrix 
and <I>+(aM) = a$+(M) for any a > and any M € M. Note that, from Caratheodory's 
theorem, a finitely-supported optimal design always exists, with m(m + l)/2 support points at 
most. We also have the following properties. 

Lemma 1 For any p G ( — l,oo), the criterion satisfies the following: 

(i) (•) is strictly concave on the set of symmetric positive definite m x m matrices; it 
is strictly isotonic on M forp € (—1,0) and strictly isotonic on M+ for p G [0,oo). 

(a) Any (f)p-optimal design ^* is nonsingular. 

(Hi) The optimal matrix M^, = M=i,[p] is unique. 



Part (i) is proved in iPukelsheimI (|l993l . Chap. 6) and the proof of {ii) is given in Appendix; 
{Hi) is a direct consequence of (i) and (ii): since an optimal design matrix M* is nonsingular, 
the strict concavity of $p (•) at M* implies that M^, is unique. Note that this does not imply 
that the optimal design measure ^* maximizing (/>p(^) is unique. 

We shall only consider values of p in (— l,oo) and, from Lemma[l]-(M), we can thus restrict 
our attention to matrices M in M+. ^p{-) is differentiable at any M G M^, with gradient 

V$+(M) = 1 [cI>+(M)]^'+iM-(^+i) = ^^M-(^+i) . 

The directional derivative ¥^^{5,;^) = \\m.^^Q+{l/ a){4>p[{\ — a)^ + av] — 4>p{Cj} is well defined 
and finite for any ^ € and any € H, with 



v) = tr{[M(z.) - M(e)]V$+[M(0]} = 



/^.xTM-(P+i)(e)xi/(dx) 1 



We shall denote by i<'^p(^,x) = F(j,^{^; 5y^) the directional derivative of 0p(-) at ^ in the direction 
of the delta measure at x. 

The following theorem, which relies on the concavity an d di fferentiability of (Oi is a 



classical result in optimal design theory, see, e.g.. iKiefed (| 19741 ) and 



PukelsheimI (jl99i Chap. 7) 



Theorem 1 (Equivalence Theorem) For any p G (— l,oo), the following statements are 
equivalent: 
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(i) ^* is (pp- optimal. 

(ii) x^M-(P+i)(^px < tr[M-P(^;)] for all x G JT". 
(Hi) ^* minimizes maxxg,^);- x) with respect to & H+. 

3 A necessary condition for support points of 0p-optimal designs 

3.1 A lower bound on x^M~'^^+^)x for the support points of an optimal design 

Take any p € (— l,oo) and any ^ G H"^. We shall omit the dependence in ^ when there is no 
ambiguity and simply write M = M(^), (pp = 4'piO- We shall also denote 

t = tiC,p) = tr[M-?'] , t^ = t^{p) = tr(M;P) , 

with M* the optimal matrix satisfying (p^ = $+(M=i,) = maxj^gs $+[M(i/)]. Define 

e = e{i,p) = niax{x^M-(P+^)x} -t. (3) 

The concavity of $+(•) implies that 0p <</>*< (/"p + F^^{^;^*) < (l)p{Cj + max^.^^- F^^{£,,yi) , 
with ^* denoting a (/)p-optimal design measure; that is, 

(^P<<p;< (Pp (1 + e/t) , (4) 

see 

Since x^M^^^+^^x < t + e for all x € see we have 

tr[M*M-(P+i)] < t + e. (5) 

On the other hand, the optimality of ^* implies (see Th. [l]-(ii)) 

tr[MM;(P+')] < t, . (6) 

Moreover, any support point x^, o f £* sa tisfies xJM* ^^^^^x* = t*. We use a construction similar 
to that in Harman and Pronzatol ((20071) and define H = H(^,p) = M-(p+i)/2mS+^M-(p+i)/2. 



Then we can write 

x7m-(p+i)x, = x7M"(p+i)/2h-V2hh-i/2m-(p+i)/2x, > Ai x7m;(p+^)x, = Ai , 

with Ai = \i{i,ip,p) = Amm(H), the minimum eigenvalue of H. Notice that Ai > 0. Ai depends 
on which is unknown. Below we shall construct a lower bound Ai on Ai and thus obtain a 
necessary condition for support points x* of in the form: 

xTm-(p+i)x, > Alt,. (7) 
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When p = (D-optimal design), we have t = t^, = m, and this necessary condition is simply 

x^M-^x, > Ai m (p = 0) ; (8) 



it corresponds to the case treated in iHarman and Pronzatd (|2007l ). When p ^ 0, is usually 
unknown and we shall use 

xJm-(p+i)x* > Xit{l + e/t)-P forp>0, (9) 
xTm-(p+i)x, > Alt for-l<p<0, (10) 

see dH) and the definitions of t, t^, (j)p, (p*. Next section is devoted to the construction of the lower 
bound A]_, using the inequalities ^ and ([6]). 

3.2 Construction of the lower bound Ai 

The inequality ([5]) can be rewritten as tr(H.^^^P~^^^'M.~^) < t + e and ([6]) can be rewritten as 
tr(H^^M~P) < t*. Consider the spectral decomposition H = SAS^, with SS"*" = S^S = and 
A the diagonal matrix whose diagonal elements are the eigenvalues Aj of H sorted by increasing 
values. Denote B = S^M~'^S and bi = {B}jj its diagonal elements, i = 1, . . . ,m. B has the 
same set of eigenvalues as M"^ and 

< 6i = A„,in(M-P) < < A„,ax(M-P) , i = l,...,m, (11) 



Magnus and Neudeckeij (|l999l . p. 211). We then obtain that (0) and © are respectively 



see, e.g., 
equivalent to 

m 

Y.b^Xy^'^'^ < t + e, 

i=l 

m 

^bi/Xi < U. (12) 
1=1 

Remark 1 Inequality (jl2p implies that Ai > 6i/t* > fei/t*. When plugged in ([7]), it gives 
x^M'^P^-'^^x,,, > bi. Although this bound is rather loose for m > 2, it cannot be improved when 
m = 1. Indeed, m = 1 implies bi = bi = t and the inequality x^M^^'^-'^-'^^x* > t is the tightest 
we can obtain, see Th. U\-(ii). In the following we shall suppose that m >2. 

Denote = X^^^^^^^ for i = 1, . . . , m > 2. The determination of Ai amounts to the solution 
of the following optimization problem: minimize wi with respect to co = (cji, . . . ,0;^)^ under 
the constraints < tJi < a;2 < • • • < ujm, Yl^i ^i^i ^ t -\- e and Yll^i ^j/'^f^^ ^ This is a 
convex problem, with Lagrangian 

(m \ / m \ 
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Its stationarity with respect to lo indicates that the optimal solution satisfies cjj = u>2 for all 
i >2. Since X^I^i = tr(M~P) = t, from the Kuhn- Tucker conditions we obtain 

biUJi + {t - bi)uj2 = t + e, 

or equivalently 

QWi + (l-a)w2 = 1 + /3, (13) 

V (1 - a)K+^ = 7*, (14) 

where a = bi/t, f3 = e/t > and 7* = t* /t. 

When p = (D-optimal design), then a = l/ m, 7^, = 1 and (fTSl). (fT^ can be directly solved 



for uJi,uj2, yielding Ai = cji to be used in ([8]), see lHarman and PronzatJ (j2007l ). However, when 
p ^ 0, a depends on M* and 7* depends on and are thus usually unknown. We must then 
determine the lowest value of uji < 0J2 satisfying ()13p . (jl4p given the information available on 
a and 7*; that is, respectively, ([TT]) which gives 1 > a > hi_/t = Amin(]V[~P)/tr(M~P), and (H]) 
which implies that 7* satisfies 

7* G [(l + /3)-P,l] ifp>0, (15) 
7* e [l,(l + /3)-P] ifp<0. (16) 

The solution to this problem is given in appendix and yields the main result of the paper. 

Theorem 2 For any p € ( — 1, 00) and any design ^ £ any point x,, € ^ such that 

x7M-(f+i)(e)x, < C{i,p)=u{-^^B{t,e) (17) 

cannot he support point of a (j)p-optimal design measure where we denoted t = tr[M~P(,^)], 
e = maXxg^- x^M~(P+-^)(^)x — t, B{t,e) = t min{l, (1 + e/t)"^}, and where uji is the unique 
solution for 9 in the interval ((0/7)^/'^^+^^, (l/7)"'^/(*'+"'^-'] of the equation 

a (I- a)P+2 

Fie; a, e, t, 7, P) = ^ + iile/t-cey-^ - ^ = ° ^''^ 

with a = Amin[M-P(e)]/tr[M-P(^)] and 7 = max{l, (1 + e/t)-P}. 

In the special case when t^, = tr[M~P(^*)] is known (thus in particular if p = 0), one can 
take B{t,e) = t^ and 7 = 7^ = t*/t in ([IT]), (fTHjl . 

Remark 2 
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1. When p is integer, F(0;a,e,t,'j,p) = is a polynomial equation in 9 of degree 2(p+l). 

2. From the definition ([3]) of e = e(^,p), 5 = maxxgx -^^ipC^jx) = e(f)p{^)/t, see and (fTTl) 
is equivalent to F0p(^,x^,) < /ip[M(^), 5] = (/>j,(^) [C{^,p)/t - 1]. Note that C{£,,p) < t, so 
that all points x such that F^^{S,,x) > are potential support points of 

3. Suppose p > with unknown and e — t- oo; then, B{t,e) 0, so that C{(,,p) — > and 
the condition ()17p brings no information on the support of ^* . The same is true when 
p < with t^ unknown and e — )• oo; 7 00, so that wi — > and again C{(,,p) — > 0. 
Suppose now that is known. Then, C{^,p) = G (Amin[M~P(^)], tr[M-P(^)]] and 

^ a/7* = Amm[M-P(0]/t* as e ^ 00, see so that C{C,p) ^ 61 = Amin[M-P(0] 
and we recover the same bound as in RemarkU^ 



4- Using a construction similar to that in lHarman and Pronzatd \200'A , Th. 3), one can show 



that the bound (jl7p with B{t,e) = t^ and 7 = t*/t gives the tightest necessary condition 
for support points: for any m > 2, any e > and any S > 0, one can exhibit an example 
with a design space ^ , a design measure ^ such that maXxg,3i:-{x^M^(P+-'^)x} —t = e, and 
an optimal design ^* with support point x* such that xJM^^^^^-'x* < co^^^ t* + (5 (with M 
and M.^ diagonal and H having eigenvalues Ai < A2 = • • • = Xm)- 



4 Examples 

Example 1. Consider the linear regression model with x = x(s) = (1, s, s^)"'", s € [—1, 1] 
(m = 3). For any p £ (— l,oo), the i;^p-optimal design on [—1,1] is unique and is supported at 
the three points {—1,0, 1}. For symmetry reasons, it corresponds to 

Ct = t6-i + {1- 2t)6o + t6i 

for some particular r* = T*{p), with r*(-l/2) = 0.45, r*(0) = 1/3 (D-optimal design), r*(l) = 
1/4 (^-optimal design) and, in the limit p — )• 00, t*(oo) = 0.2 (i?-optimal design), see Fig.[2]-left 
for a plot of T*{p) for p € [—1/2, 1]. Here, 5s denotes the Dirac delta measure at s. 

To illustrate the impact of not knowing on the construction of uj^^^ through the solution 
of (fT8|) . we take p = 1 and compute w^"*"^ for the cases 7 = 1 (t* unknown) and 7 = t*/t 
(t* known) for different designs S^r, t G [t*(1) — 1/16, t*(1) + 1/16]. Figure [1] shows that the 
value obtained for t^, unknown (solid line) is not much worse, i.e., smaller, than the value for 
t* known (dashed line). Note that considering different designs with r 7^ T*{p) is equivalent 
to considering different e given by ([3]), with e being approximately linear in |r — T*{p)\ for the 
range of values of r considered. 
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Figure 1: Value of uj^^^ for different designs ^-r, t G [3/16, 5/16] (p — 1, t, unknown in solid line, t, known in 
dashed line). 

The marginal deterioration of the bound (117p due to the ignorance of when e is small 
enough is further illustrated by Fig. [2j Here, we set e at some fixed value (the values e = 0.1 
and e = 1 are considered), and for values of p in the range [—1/2,1] we compute T{p,e) such 
that max^e[_i,i] xT(s)M-(p+i)(^^)x(s) = tr[M-P(^^)] + e. The values of T*{p) and r(p,0.1) are 
shown in Fig. [2]-left, in solid and dashed lines respectively. Then, for each p and associated 
design Cr{p,e) we compute the bound C(^^(p^£),p) of (fT7|) in the two situations unknown and 
known; see the plots in Fig. [2]-right. Note that bound for t^, unknown (solid line) remains near 
the bound for known (dashed line) when e = 0.1; the situation deteriorates for larger e but the 
two bounds get close as p approaches and exactly coincide at p = (since then t = t^, = m) . 



Example 2. Take now the complete product-type interaction model x(s) = x(si) (8) x(s2) 



s = (si, S2), with (8) denoting tensor product and x(si) = (1, Sj, s?) , Sj G [—1, 1], for i = 1,2 
(m = 9). The D-optimal (respectively ^-optimal) design for this problem is the cross product 
of two D-optimal designs (resp. ^-optimal designs) for one single fact or, i.e., it corres ponds to 
the cross product of two designs ^t- with r = 1/3 (resp. r = 1/4), see ISchwabd (119961 . Chap. 4 
and 5). The optimal values of 4>p, p = 0,l, are (/>S = 16^/^/9 ~ 0.2800 and (pl = 9/64 ~ 0.1406. 
We consider the iterative construction of optimal designs through the recursion 

[xTM-(P+i)(efc)x,]» 



fc+1 



Ef=\[x7M-(^'+i)(e,)x,]« 



(19) 



where the design measure at iteration k, allocates mass wf at the point Xj present in ^ at 
iteration k, i = 1, . . . ,Nk. The initial design space corresponds to a uniform grid for s, with Sj 
varying from —1 to 1 by steps of 0.05 (41 values), i = 1,2, which gives A'o = 1681. The initial 
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Figure 2: Left: T*{p) such that Ct*(p) ~ is </>p-optimal for p {solid line) and r(p, e) such that 
max3g[_i 1] x^(s)M~'*'+-''(5^(p_j))x(s) = tr[M~''(^^(p ,,))] + e (e = 0.1, dashed line). Right: bound C(^T-(p_j),p) in 
(|17[1 for the two cases t* unknown {solid line) and known {dashed line) for e = 0.1 and e = 1. 



design is the uniform measure on those 1681 points. We take a = 1 for D-optimal design 
(p = 0) ^-iid g = 1/2 for A-qptiin al de sign (p = 1), wh ich ensures monotonia con vergence to the 
optimum, see iTitteringtonI (jl976l ) and iPazmanI (jl986l ) for D-optimal design and iTorsnevI (|l983l ) 
for A-optimal design; see also Fig. [HJ-left. Due to the convergence of to the optimal design, 
Cfc = ^{^k) given by ^ is decreasing with k, see Fig. [31-right. 

We use inequality (fT7|) to reduce the cardinality of ^ when possible: any point that 
violates (|17p cannot be a support point of the optimal measure and is removed from ^ . Here we 
simply set its mass to zero and rescale the weights of remai ning point so that they sum t o one, 
but more sophisticated reallocation rules can be used, see iHarman and Pronzatd (|2007l ) . Nk 
thus decreases with k, rendering the iterations (jl9p simpler and simpler as k increases. Figure H] 
shows the evolution of with k, both for D-optimal and A-optimal designs. The decrease of 
Nk is slower for the latter, the bound C{S,,p) in (fT7|) being more pessimistic, see Fig. [2l-right, 
and e being larger, see Figl3]-right. Note that the cancelation of points does not hamper the 
convergence of (jl9p since ()17p is used a finite number of times only (obviously bounded by A'^o) 
— the heuristic rule used to reallocate weights of points that are removed may, however, impact 
monotonicity, although this is not the case in the present example, see Fig. [3]-left. 



5 Possible extensions and conclusions 

Multivariate regression and Bayesian optimal design involve information matrices that can be 
expressed as M(^) = J^. .'#(x) ^(dx) with ^(x) G M having rank larger than one (we sup- 
pose that ^(•) is measurable and that {^(x),x € ^} forms a compact subset of M). The 
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Figure 3: 4'i£,k) — left — and = e(Cfc) given by ^ — right — as functions of k for the recursion (|19p : 
D-optimal design is in dashed line, A-optimal design is in solid line. 



1800 




10° 10^ 10^ lo'' 

k 



Figure 4: A'^^ as a function of k when using (|17[l to remove points from J?^'; for _D-optimal design {dashed line) 
and ^-optimal design {soUd line). 
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results presented here can ea sily be extended to that situation, following the same lines as in 



Harman and Trnovskal (120091 ) where the case p = is considered. 

The i?-optimality criterion (pEiO = = Amin[]V[(^)] is not differentiable in general, 

but is differentiable at M when Amm(M) has multiplicity one, with gradient Vi?!>£;(M) = 

vv"'' where v denotes the eigenvector of unit length (unique up to a sign change) associated with 
Amin(M). Although 4>E{i) Corresponds to the limit of 0p (C) as p tends to infinity, the results 
of Sect. [3] do not extend to this limiting situation, even in the differentiable case; i?-optimality 
thus requires a special treatment and will be considered elsewhere. 



The determination of a - D-optim a. 



desi gn ca n be u s ed for maximum-likelihood estimation 



in mixture models, see, e.g., iLindsayl (119831 ) and iMalleti (119861). and for the construction of the 



Khachivan and Todd 



minim um- v olume ellips o id con taining a compact set, see, e.g.. lSibsonl (|1972| ). 
(|l993l ) and iKhachivanI (jl996l ). More generally, for any q G (— l,oo) the determination of the 
eUipsoid (f (A) = {z G M'" : z'^Az < 1}, A G M, containing the k points xi, . . . ,Xfc of 
and such that 0g(A) is maximum is equivalent to the determination of a i?i>p-optimal design on 
^ = |xi, . . . , Xfcl wit h p = —q/{l+q) G (—1, oo), and the optimal matrix A.,, equals M^. ^^"""^^/t*; 
see iPukelsheimI (|l993l . Chap. 6). The delimitation of the support points of a (/>p-optimal design 
can therefore also be used to accelerate the algorithmic construction of "(/)|j-optimal ellipsoids" 
containing compact sets. 

In Sect, m we considered the suppression of points that cannot be support points of an 
optimal design in a multiplicative algorithm. When J^T is not finite, or is finite but very l arge, it 



Wu 


(1978 


) and 


Bohnine 


(1986) 



fe, of a point Xfc of ^ that maximizes F(^p(^fc,x) given by ([2]), at least approximately. Condition 
(|17|) of Theorem [2] can then be used to restrict the search for a suitable x^, in a domain that 
shrinks as k increases. Further developments are required to construct algorithms making an 
efficient use of ()17p for the inclusion of new support points. 



Appendix 

Proof of Lemma{l\-(ii). 

For p >0 the result follows from the observation that (M) = when M is singular while 
there exists a nonsingular 1V[(^) with [M(^)] > 0. 

For p G (—1,0) we prove the result by contradiction. Take any Mq = M(,^o) singular and 
consider its spectral decomposition in an orthonormal basis of eigenvectors Vj, i = 1, . . . ,m: 
Mq = Yl^i ^i"^! ■ Suppose that the eigenvalues Aj are sorted by increasing values, so that 
Aj = for i = 1, . . . , s when the eigenvalue has multiplicity s (and Mq has rank m — s). Since 
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there exists a nonsingular design, ^ spans M™ and each eigenvector Vj, i = 1, . . . , s, can be 
written as Vj = 7j /^u(_5r) z/ij(dz) for some 7^ > 0, where ^j(-) is a probabihty measure on 
the compact set U {—^) (in fact, from Caratheodory's theorem, one may consider finitely 
supported measures only, with m + 1 support points at most). Then, 



7i 



zz^//i(dz) - Viv] 



(7i z - Vi)(7j z - Vj)"^ Aij(dz) , 



which is non-negative definite. Denote ^(dz) = 7^ Mi(dz)]/(^^^j^ 7?), which defines a 

probability measure on ^ U {—^). We thus obtain that, for any a G (0, 1), the matrix 



(1 - a) Mo + a 



I3 



zz^ li{dz) 



(1 - a) Mo + a 



E 



ViV, 



is non-negative definite. Now, f^^^_^~^zz~^ fi{dz) can be written as xx^ /i(dx), where 
fl{s^) = fi{£/)+fj,{—£/) for any measurable set £/ C and /!(•) is thus a design measure on 
Therefore, $+(Ma) > (M^) , where M^ = M[(l-a)^o + a/2] = (l-a)Mo+a xx"^ /I(dx) 
and M^ = (1 — a) Mq + (a/p) X]i=i^*^j'^' with p = Yli=i^i- "^^^ eigenvector decomposition 
of Mq gives 



*;(ML)= - 



(l-o)-" 




+ a 



-i/p 



which reaches its maximum value for a = a* = [1 + {p ^ Ei^^+i ''/■s)"'^^*-^^^-*] ^ G (0, 1). It 
implies that $+(Mq*) > $+(M^.) > $+(Mo); that is, 4>p[{l - a*)^o + a*fi] > 0p(^o) and ^0 is 
not optimal. ■ 

Proof of Th [H The proof is in three parts. In (i) we show that for given a and 7* the 
equations p^ . with uji < UJ2 have a unique solution a;J(a,7*) for wi, with a;|(a,7*) G 
((a/7*)"'^/^^'''^^, (1/7=|,)^/*-^'''^^]. Then in (iz) we show that this solution is non-decreasing in a, so 
that the required lowest bound is obtained for a = bi/t, see pT|) . Finally, in (iii) we consider 
the case when is unknown. 

(i) Expressing UJ2 as a function of coi using (fT3l) . we obtain 0J2 = /i(<^i) = (1+/3— awi)/(l— a), 
i.e., a decreasing linear function of oji with slope — a/(l — a) and such that + /3)/a] = 
0. Doing the same with (jl4p . we obtain a;2 = f2{^i) with /2(-) decreasing and concave for 
UJi G ((a/7*)^/*^^^^\ 00), /2(0) tending to infinity when 6 approaches (0/7*)^/*^^+^^ from above 
and lime_,oo/2(^) = 1/a - 1- Note that (USD, ^ imply that {a/-f,y/^P+^^ < (1/7*)^/(p+^) < 
(1 + /3)/a. Therefore, f2{9) > fi{G) for close enough to (0/7*)^/'^^'*'^) or large enough. 

Denote f2{0) = df2{9)/d6 and consider = (1/7*)^/'^^+^). Direct calculations indicate that 
72(6'*) = 0*, f2{^*) = — Oi) with, moreover, fi{9^) > f2{9*) when /3 > 0, i.e., when e > 0, 

due to US]) and Two solutions a;| ^, t^jf thus exist for ([T3D, ([H]), with wjf ^ < 6** < tj* ^. Only 
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^ is such that the associated ^2 a satisfies a > a- When e = 0, then fi{0*) = f2{d*) = 0* 
and the two solutions a;*^,a;^j, are confounded and equal 0* (and also coincide with uj2a and 
uj2b)- '^^^ equations (fT3]) and with wi < UJ2 thus always have a unique solution a;|(a,7*) 
and this solution belongs to the interval ((a/7*)"^/*^^'^"^\ 0*]. 

(ii) Applying the implicit function theorem to p^ . (fT^ we obtain that the solution wf (a, 7^,) 
satisfies 

a(p + l)[(^)P+2 - (w*)P+2] 



[{p+l){l-z) + z{zP+'-l)], 



a(p + l)(zP+2 - 1) 

where 2 = ^2/^1 > 1- Denote /(z) = — — its derivative is d/(z)/dz = {p+ 

2)(zP+i-l) so that /(z) > /(I) = 0. Since (HU gives a > 61/*, one hasa;^(a,7*) > 7*)- 
The substitution of 7*)]^+-*^ for Ai in ([7]) concludes the proof for the case when is 

known. 

{in) When is unknown, an upper bound can be substituted for in (fT2]) . Using (fT5]l . 
p^ . this amounts at replacing 7^, by the upper bound 7 = max{l, (1 + e/t)^^}. The necessary 
conditions Q, ([TO]) with Ai = [a;f (61/*, 7)]^"^^ then give ([E 
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